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Abstract. 

We introduce a model for yielding, inspired by fracture models and the failure of 
a sheared granular medium in which the applied shear is resisted by self-organized 
force chains. The force chains in the granular medium (GM) are considered as a 
bundle of fibres of finite strength amongst which stress is randomly redistributed 
after any other fibre breaks under excessive load. The model provides an exponential 
distribution of the internal stress and a log-normal shaped distribution of failure stress, 
in agreement with experimental observations. The model displays critical behaviour 
which approaches mean field as the number of random neighbours k becomes large and 
also displays a failure strength which remains finite in the limit of infinite size. From 
comparison with different models it is argued that this is an effect of uncorrelation. 
All these macroscopic properties appear statistically stable with respect to the choice 
of the chains' initial strength distribution. The investigated model is relevant for all 
systems in which some generic external load or pressure is borne by a number of units, 
independent of one another except when failure of a unit causes load transfer to some 
random choice of neighbouring units. 



PACS numbers: 62.20.M-,64.60.av,45.70.-n,05.65.+b 



Submitted to: Journal of Statistics and Mechanics: theory and experiment 



Random neighbour model for yielding 



2 



1. Introduction 

Many different processes, from stress propagation in a heterogeneous medium to the 
services for client satisfaction involve redistribution of loads when some units break 
down, potentially triggering an avalanche of failures and causing catastrophic breakdown 
of the entire system. Other examples might be data centers which supply information 
over some communications infrastructure and electrical power supply from power 
stations. Such generic non-linear and irreversible failure processes, where the internal 
dynamics can be highly susceptible to external perturbations, constitute a general 
feature of yielding processes and may be observed in many other different contexts and 
systems such as fracture, plastic deformation, granular slip, economic crashes, etc. This 
class of phenomena generally exhibits critical behaviour with relevant characterizing 
features such as finite yield points and broad avalanche amplitude distributions, as well 
as a highly inhomogeneous internal load distribution over the load-bearing units. 

In this article, we propose a simple model for yielding that is general in purpose but 
essentially inspired by the case of a Granular Medium (GM) under shear stress. The 
distribution of stress within a GM is a characteristic feature at variance with fluids; it 
is neither isotropic nor homogeneous due to the presence of force chains [U [2] , highly 
singular random paths which bear most of the load and screen the majority of grains from 
stress. In real GMs, the failure of a force chain triggers stress transfer processes which 
involve complex spatial structures [3] which can also dynamically change under the chain 
rupture and alter the redistribution process itself. These cannot easily be incorporated 
in any simple model without resorting to molecular dynamics simulations. The model 
presented here uses a simplifying approach in which force chains, each represented by 
a single number, are the smallest components. Each chain in the medium is able to 
support a maximum load, and as the increasing external stress causes the weakest 
remaining fibre to break, its load is redistributed to a number of other surviving chains. 
The model shares some features with the Fibre Bundle Model (FBM) for fractures [I], 
but introduces a novel random load redistribution mechanism. For this reason, it shall 
henceforth be denoted the Random Fiber Bundle Model (RFBM). With respect to 
the sheared granular bed, fibres play the role of the force chains within the granular 
material, an analogy already exploited in [5] for compressional stress in which broken 
chains undergo a restructuring which results in a stiffening of the bundle. In the RFBM 
model presented here, representing brittle shear failure, breaking force chains transfer 
their stress to some random selection of force chains. The model aims to reproduce some 
macroscopic statistical features of the yield process already known from experimental 
measurements, as follows. 

i) Experimental studies [SI U\ have pointed out that stress distributions within 
a granular medium are predominantly characterized by long exponential tails. This 
property seems to be general since it appears in different regimes of granular mechanics 
ranging from static assembly [6] to dynamic shear [8j. A simple explanation for the 
static case is given by the so called q-model [9] which, although imprecise in some 
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aspects [ID] , produces robust exponential tails. In the q-model, grains arranged in 
horizontal layers discharge random proportions of their load on grains in the layer below. 
The model and its developments have been considered and discussed in several works 
[TT| [T2T ITDl IT3T EH EH], from which the idea emerges that random load transfer can 
produce long exponential tails as a stable characteristic. The RFBM model proposed 
here allows force chains to irreversibly break, but retains this random rule for load 
redistribution, by randomly selecting which of the remaining fibres should bear the 
load of the breaking fibre. The randomness of the choice may be interpreted as an 
expression of our inherent ignorance about the real stress and complex spatial force 
chain correlations. 

ii) Another characteristic of the shear yield process is that the force required to 
trigger a slip in a sheared granular material, as well the force required to sustain steady 
motion, are statistical quantities that have been shown to accurately follow log-normal 
shaped distributions [161 EI]- This feature does not seem to be exclusive to granular 
yield and has in fact been observed in solid-on-solid [18J and fibre-on-fibre [19] friction 
not to mention more generic yielding phenomena ranging from friction and wear [20] 
to cohesion and bearing capacity of soils [2T] and the strength of wood structures [22J. 
Moreover, log-normal shaped distributions may also be produced by a variety of model 
systems ranging from fracture [23] to dislocation [23], as well as magnetic and self 
organized systems [25J. 

Hi) Finally, many yield processes display critical features. Typical examples are 
the depinning transition of elastic manifolds [26] and fracture processes [26]. Granular 
dynamics proper has been shown to display critical dynamics in the chaotic stick-slip 
motion phase [27] . 

It will be shown that the RFBM reproduces all the main features described above, 
i.e. i) the exponential distribution of internal loads, ii) a non Gaussian (log-normal 
shaped) probability distribution of failure loads and Hi) a critical dynamics characterized 
by a mean field exponent and a cut-off, the divergence with size of which depends on 
the number of abssorbing chains. In addition, it will be seen that the model exhibits 
a finite yield strength even in the infinite system limit. Although this latter case is 
difficult (sic!) to test experimentally, it is at least clear that it holds in finite size cases 
[TBI |28| 129] . From the comparison with models having different geometries and load 
sharing rules we will conclude that such a feature is related to the absence of spatial 
correlation due to the random neighbour allocation. 

This article is organized as follows: Section [2] describes the RFBM and recalls some 
related models in order to compare their properties. The distribution of the internal 
stress is investigated in section [3] and section H] discusses the statistical properties of the 
yield point stress. Section [5] shows the critical aspect of the process, drawing a careful 
comparison with those of the usual FBM. Concluding remarks are drawn in section [63 
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2. Model details and features 

The dynamics of the model discussed in the present work is to some extent based on 
the Fibre-Bundle Model (FBM), a simple but revealing approach to the study of failure 
originally introduced by Peirce [4] to explore the properties of textiles and subsequently 
extensively explored [3 EQl ED E21 E31 ES EH]. The FBM consists of a bundle of N 
harmonic fibres loaded in parallel by a slowly increasing external load F. Each fibre 
has a different breaking strength Sj drawn from some identical initial distribution t(s). 
As the external load is increased, fibres whose load exceeds their strength break, and 
redistribute their load equally on each of the remaining intact fibres R. ("global load 
sharing" law). Every intact fibre therefore, carries the same load and in this respect it 
can be considered the simplest mean field model for fracture. 

As in the FBM we define the RFBM by a set of harmonic fibres each having 
identical elastic modulus but being able to support a different random maximum failure 
load Si extracted from some probability distribution t(s). The bundle is subjected to an 
external load F, initially zero, which is then increased so as to break only the weakest 
fibre which, in the RFBM, then redistributes its load on a fixed number k of different 
randomly chosen intact fibres, any of which may now exceed their capacity and break 
in turn. Thus we may observe an "avalanche" of breaking fibres. After the avalanche (if 
any) the external load is increased again so as to break only the new weakest fibre (this 
corresponds to the limit of adiabatic forcing). We investigate several variants of this 
basic RFBM. The case k — 1, in which a breaking fibre passes its entire load to a single 
other fibre; the case k = 2 where the load is redistributed in i) equal proportions (ES) 
or ii) randomly chosen proportions p and q = 1 — p (PQS) among the two randomly 
chosen fibres (p is newly extracted at the breaking of each fibre). For the cases k = 3, 4 
some main properties are discussed. The rules of the RFBM bear some similarities with 
the conservative limit of the random neighbour sand pile model [361 [37] except that here 
the threshold of each fibre (site) is drawn from a random distribution and fibre breaking 
is irreversible. However, like the sandpile model, every intact fibre bears a different load 
which depends on the fibre's history. 

We define the stress as a = ^Ti/N = F/N, where Tj is the load on the i^ 1 
fibre and N is the initial number of fibres. The yield stress a c then, is a evaluated at 
the beginning of the final avalanche which ruptures the bundle as a whole. We define 
the reduced bundle load as / = ojo c which varies from zero to one as the simulation 
proceeds, and permits the comparison of results from different realizations. 

In the original FBM, a c is a random number for small N, with a skewed distribution 
which tends to a Gaussian with increasing N and a 5 function for iV = oo. A non- 
mean field version of the FBM also exists, with fibres placed on a regular lattice, and 
discharging to their nearest neighbours only (Local Sharing rule - LFBM) [381 EI]- The 
behaviour of this model has also been recently investigated on different types of complex 
networks (CFBM) [39]. For comparison we will also present and discuss some features 
and results from these variations. It will be seen that the RFBM exhibits behaviour 
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which in some sense places it between the FBM and LFBM. Except where otherwise 
specified, all simulations are performed on systems with N = 10, 000 fibres. 

3. The internal distribution of loads 

As noted in the introduction, our main motivation for devising and investigating the 
RFBM is to attempt a description of the yield of a sheared granular bed. In doing so, one 
reduces the granular force chains to simple fibres, characterized by a single variable s, 
their ability to bear stress. It is known that a GM under either static or dynamic loading 
displays an exponential distribution of internal forces [6j El [TTJ [7J E2l [TOj [131 EH ESI E] , 
a property generally encountered in disordered systems [10, SH El] ■ We have computed 
therefore the internal load distribution p(r) as a function of the reduced bundle load 

/ = (?/ Ve- 
in figure [U we demonstrate this distribution at the yield point / = 1 for four different 
choices of the fibres failure threshold distribution t(s) (uniform, exponential, Gumbel 
and Weibull, details in the appendix), with the random neighbour number k = 2 and 
equal sharing (ES). The figure shows that the RFBM exhibits robust exponential tails 
in the load distribution. Since initially / = and the load on each fibre is zero, this 
distribution builds up according as the simulation proceeds and the external load is 
increased. Figure [2] illustrates the evolution of the internal load distribution p{r) for 
increasing values of / and t(s) exponential. The ES case is shown in a), whereas b) 
shows results for random sharing on k = 2 fibres (PQS). It is clear that the exponential 
tail in the distribution of internal forces is achieved quite early in the evolution of the 
system, certainly long before yielding occurs. All variants of the model investigated 
display exponential tails similar to the k = 2, in particular for k = 1,3,4, and using 
uniform t(s) instead of exponential. 

From the above cases we conclude that exponential tails in the internal distribution 
of loads are robust features of the model, and the generality of results when changing t(s) 
and k suggests that the random load redistribution property of the model is a general 
mechanism naturally leading to the exponential load distribution [13], analogous to that 
of the g-model. As a test we have also computed p(r) for the LFBM in a ID periodic 
geometry, in which the load of a broken fibre is transferred to k = 2 neighbouring fibres. 
The resulting p{r) at different values of the external load are shown in figure [31 In both 
cases investigated (t(s) uniform and exponential), the internal load distribution displays 
long algebraic tails. This, in itself, is a remarkable feature of the LFBM and merits 
further study. From these results and those related to the q-model and similar models 
[13] . we conlcude that the exponential distribution of loads obtained above is specific 
to random correlations. This, of course, does not exclude the possibility that other 
mechanisms may also generate exponential distributions of loads, though we suggest 
that any other such mechanism will hardly be simpler that that presented here. 
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4. The yield point 

In this section we investigate the behaviour of the RFBM at the system's global yield 
point, varying a number of parameters such as the bundle size N and the fibre strength 
distribution t(s). We compute in particular g(cr c ), the probability that the system yields 
at external stress a c . Among other results we find that: i) the dependence of the shape of 
g(cr c ) on the fibre strength distribution t(s) is very weak; ii) for not too large N, g{o~ c ) 
has a log- normal shape and Hi) the mean, variance and skewness of g(o~ c ) decreases 
algebraically with the system size. 

We consider initially the case k = 2. Experimental measurements of the yield point 
of a sheared granular medium have been performed extensively in [T6l [T7] . and have 
demonstrated the distribution to be skewed and far from Gaussian. It seems remarkable 
that it is also well fitted by a log- normal distribution (or similar) [16] . Such a distribution 
seems to be a generic signature of the strong correlations introduced into the system by 
the force chains [TT] , and some authors have resorted to universality as some underlying 
feature for such observations [25J. The RFBM model outlined here produces a yield 
stress distribution g(a c ) wich is also very well approximated by a log- normal: 

2 



9( a c) = — — : — exp 

V27r \a c — 0o S\ n 



ln(cr c - (T ) - AUn 



y/25 



in 



(1) 



where /ii n and 5\ n are the mean and standard deviation of ln(cr c — <7o) an d o"o is a rigidity 
threshold which imposes a lower bound to the minimum yield stress obtainable in any 
finite system realization. 

Figure H] shows the pdf g(c c ), for N = 100 fibres and four different choices of the 
fibres' initial strength distribution t(s): uniform, exponential, Gumbell and Weibull 
(see the appendix for details on the precise forms of t(s)). Both ES (continuous lines) 
and PQS (dashed lines) are shown, each result averaging over 10 6 realisations. When 
a standardized variable is adopted z = (ln(a c — a ) — /i ln ) /<5 ln all curves collapse to the 
same normal curve, showing that g(a c ) is not strongly affected by the choice of fibre 
strength distribution t(s). For very low N, therefore, g(o~ c ) is far from Gaussian [30J, 
becoming Gaussian only for sufficiently large N. This also follows from the expression 
for z since: /ii n = (ln(cr c — <7o)) = ln(/i) where jl = \/Y[( a c — o"o) is the geometric mean 
of (cx c — o"o) over n realisations. As is the case with large N, the distributions is narrow, 
and the geometric mean approximates the arithmetic mean jl — > (a c — cr ) = (a c ) — a 
and so 

* = M*c - *„) - hi«a c > - a Q )} = f In ( 1 + ^) (2) 

Thus, if < 1 (i.e. the distribution is narrow) then g(a c ) approximates a Gaussian 

centered at (a c ) with standard deviation A = [(a c ) — cr ] 5\ Q : 



V2A 



(3) 
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Table 1. Critical exponents for the parameters of the distribution of failure stresses 
g(cr c ), as defined in (TJ|. 



Model 


Coo 




X. 




FBM 


0.25 


-0.66 


-0.99 


-1.91 


RFBM ES 


0.20 


-0.56 


-0.99 


-1.66 


RFBM PQS 


0.17 


-0.5 


-0.89 


-1.30 



We study now the size dependence of the mean fi ac , the standard deviation 5 ac and 
the skewness r y CTc of the yield point distribution g(o~ c ), for sizes N from 10 2 to 10 6 as 
shown in figure The ES and PQS cases are shown and compared with the FBM case. 
We focus on the case of a uniform t(s) which is more commonly adopted for the FBM 
and observe that good power-law fits may be obtained: 

S Uc oc N~ x (4) 

la c OC N-* 

where indicates a non-zero failure load per fibre for an infinite bundle N = oo. 

The scaling of \i ac in figure [5] a) shows that an algebraic scaling is common to the 
FBM (global sharing) and retains a finite strength as iV — > oo. This is in constrast 
to the LFBM where, for large N, one obtains instead fi Uc ~ 1/ ln(iV) [23J, i.e. a vanishing 
strength limit for N — >■ oo. Despite the local sharing rule in CFBM models, however, 
a finite strength limit is again retained [39] as in the RFBM; this may be due to the 
random network topology that destroys spatial correlations. Thus the requirement for 
non-zero appears to be an absence of spatial correlations, due either to annealed 
randomness (RFBM) or quenched randomness (CFBM). Notably, all variants of the 
RFBM studied possess a systematically weaker rupture threshold than the FBM, as 
they tend to concentrate the external load on some subset of fibres. 

Figures [5] b) and c) show the scaling of the standard deviation and skewness of 
g(cr c ) f° r sizes N from 10 2 to 10 6 fibres. The skewness vanishes for N — > oo and the 
distribution tends to a Gaussian, consistent with the FBM (also shown) though the 
RFBM (both ES and PQS), for increasing N, gives rise to symmetric g(cr c ) somewhat 
later than the FBM model. The standard deviation S ac — > in the limit iV — > oo, as 
expected, reducing the distribution g(a c ) to a delta-function at a^. Scaling exponents 
appearing in (Referenceseq:phichipsi) for the models investigated are recorded in table [TJ 

The dependence of the mean, standard deviation and skewness of the yield stress 
distribution g(o~ c ) on the number of random neighbours k is shown in figure [6j Values 
are reported for two different choices of the system size, N = 100 and iV = 10000. They 
are shown to rapidly approach their mean-field values as k increases. 
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5. Critical features 

One remarkable feature of the FBM is that, during loading, fibres break in bursts giving 
rise to a highly intermittent and erratic activity with no characteristic size. This kind 
of activity is actually observed in real fracture experiments (e.g. through the emission 
of elastic waves), from microfractures and crystal defects to earthquakes, and is one of 
the reasons that this model is the subject of many investigations [23]. Considering all 
the bursts occurring from / = to / = 1, the probability of obtaining a burst in which 
n fibres break (i.e. 'size' n) is [42J: 

p(n) ~ n~ 5 / 2 

When looking at the same distribution for a given interval of load / ± 5f one obtains 
p(n) ~ n~ 3 ^ 2 exp(— n/n co ) (5) 

with 

n co « 1/(1 - /T; 7 = 1- (6) 
That is, a characteristic local scale exists when the system is far from the global rupture, 
which however diverges algebraically as the final failure is approached. This behaviour 
is shown for the original, mean field, FBM in figure [7| a), where p(n) and n co (inset) 
are fitted according to (jSJ) and (El), and indicates that the whole fracture process can 
be viewed as an approach to a critical point, where lim/_>.i ri co = oo with a divergence 
exponent 7 = 1, and avalanches of any size are possible. Similar behaviour is observed 
in the CFBM [39] but not in the LFBM, where the distribution observed is stretched 
exponential [381 EI]- Critical behaviour in two dimension is only observed for when 
sharing is long range enough [5] . 

In figure [7] b) the evolution of the avalanche distribution for the ES RFBM with 
k = 2 is demonstrated. As in the FBM, the initial avalanches are of course small but the 
distribution broadens as / increases. Curves are successfully fitted by (jSJ), indicating 
behaviour identical to the FBM. The inset displays the divergence of the cut-off n co 
as / — > 1 together with a fit with ([6]), again identical to the FBM model, but with a 
different exponent: 7 = 1.28. The same critical behaviour is observed for the RFBM 
with different values of k. For instance, 7 = 1.08 for k — 4. In fact, the exponent 
ruling the cut-off divergence is seen to approach that of the FBM as k increases, as 
expected. This is illustrated by figure El where the evolution of the parameter n co is 
shown as a function of 1 - / for the FBM, the ES RFBM with k = 2,3,4 and the 
PQS RFBM (with k = 2). In this last case 7 has a very high value 7 = 1.5. For ES 
the system critical behaviour rapidly approaches that of the FBM [42] and CFBM [39] 
as k increases. Critical behaviour is not observed in the LFBM with dimension = 1; 
on the other hand, CFBM and the model presented here utilise a quenched and an 
annealed random neighbour selection respectively, and so may be considered to have a 
high dimensionality. This agrees with the observation that even LFBM in dimensions 
> 2 also show criticality [5]. 
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6. Summary and discussion 

We have introduced a model for yielding (RFBM) specifically inspired by the yield 
of a sheared granular bed under increasing stress. In the model, when a force chain 
fails it discharges its load on k other randomly chosen chains. Each force chain can 
withstand a finite amount of stress extracted from some probability distribution. In this 
respect the model is closely related to the Fibre Bundle Model (FBM) but employing a 
novel, random load-sharing rule. This model may be considered a schematization of any 
general system in which some external load or pressure is borne by a number of units, 
each independent of the others except when failure triggers some load-transfer process. 

For a finite number of units N, the yield strength at which global failure occurs 
displays a probability distribution with a log-normal shape, a feature common to many 
systems. This property is shown to be largely independent of the choice of units' 
capacities distribution t(s). For increasing N the distribution exhibits scaling and 
finally shrinks to a non-zero deterministic value, at variance with the LFBM (local 
sharing rule). This implies a residual strength for the systems in the limit iV — > oo, a 
property that we attribute to the lack of spatial correlations. 

Under increasing stress and soon after the loading starts, the distribution of internal 
loads evolves to an exponential, a feature characteristic of granular media and that also 
is found to be very robust with respect to the choice of the capacities distribution 
t(s). In addition it is also essentially independent of the random redistribution rule for 
broken units. Considering then the evident analogies with the g-model, we hypothesize 
that random neighbour selection is a condition for the emergence of the exponential 
load distribution. This is in stark contrast to the power-law internal load distribution 
obtained from the ID LFBM. Considering also the results relating to random networks, 
we also conclude that the absence of spatial correlations gives rise to the critical 
fluctuations that the model exhibits as the global breakdown is approached. 
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Appendix 

List of probability densities for the threshold distribustions employed in the examples 
and related parameters. 

1) uniform: 

t(s) = 1 for s G [0, 1), t(s) = otherwise; 

2) exponential: 

t(s) = e~ s for s > 0; 
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3) Gumbel: 

t(s) = e -( s+e ~ s ) for s > 0; 

4) Weibull: 

t(s) = 2se~ s2 for s > 0. 
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Figure 1. The distribution of the load per fibre p(r) for the RFBM ES (k = 2) at 
the yield point / = 1. p(r) shows stable long exponential tails for several different 
probability distributions of the fibre strength t(s) (see Appendix for the distributions 
details). 




Figure 2. Distribution of the load per fibre p(r) in the, k = 2, RFBM for an 
exponential distribution of fibre strenght: t(s) — exp(-s), s > at increasing values 
of /: a) Equal sharing (ES); b) Random sharing (PQS). Exponential tails in p(r) are 
established early in the system's evolution. 
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Figure 4. Distribution of the yield strength z in the standardized variable (see text) for 
different fibre strength distributions t(s): uniform, exponential, Gumbel and Weibull 
(see Appendix). The distribution g(z) appears robust with respect to the choice 
of t(s). Continuous and dashed lines refer to ES and PQS respectively. The inset 
displays the same results with non-standardized variables; PQS is seen to be slightly 
but consistently weaker than ES. 
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Figure 5. Variation of the failure distributions parameters with the random neighbour 
number k for N = 100 (black) or N = 10000 (red) in the ES case, a) The mean b) 
the standard deviation and c) the skewness of g(a c ). As the number of neighbours 
increases, these parameters rapidly approach the mean field values. 
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Figure 6. The mean /v c , standard deviation 5 ac and skewness 7 CTc of the failure stress 
distribution g((J c ) are shown here as a function of the bundle size N for the FBM and 
the RFBM ES & PQS: a) average, b) variance, c) skewness. There is a clear power-law 
dependence on N for each of the variables. The skewness for N = 10 6 is omitted as it 
was not possible to obtain a reliable estimation. The exponents are reported in table [T] 
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Figure 7. Avalanche size distribution p(n) at several values of / for: a) The FBM 
(global load sharing, t(s) uniform) at several values of the external load /, together 
with the curve fits by ([3]). b) he RFBM ES (fc = 2, uniform t(s)). Insets displays the 
dependence of the cut-off n co on / according to ([5]) . 
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Figure 8. Divergence of the avalanche cut-off n co with / according to ([6]) for the 
FBM, the RFBM ES k = 2, 3, 4, and for the RFBM PQS k = 2. For increasing k the 
critical exponent 7 approaches the mean- held FBM value of 7 = 1, whereas the PQS 
displays 7 = 1.5. 



